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Abstract 

We derive analytical expressions for the solid angle subtended by a circular disc 
at a point source with cosine angular distribution (/(//) = /V 71 ") under the sole 
condition that the disc lies in the half-space illuminated by the source (fi > 0). The 
expressions are given with reference to two alternative coordinate systems (S and 
S'), S being such that the z axis is parallel to the symmetry axis of the disc and 
S' such that the z' axis is aligned with the source direction. Sample plots of the 
expressions are presented. 
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1 Introduction 



In two recent works we obtained analytical expressions for the solid angle sub- 
tended by a cylindrical shaped detector at a point cosine source in the cases 
where the source axis is orthogonal [1] or parallel [2] to the cylinder axis of 
revolution. As ancillary results we also derived expressions for the solid angle 
denned by a circular disc in the cases where the source axis is orthogonal (f2j_) 
or parallel (On) to the symmetry axis of the disc. This latter result ap- 
peared in a previous work by Hubbell et al [3, eq. 29], where it is also credited 
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to other authors [4,5]. In that work, a quite general treatment of the radiation 
field due to a circular disc source with axial symmetry is given in terms of 
a Legendre expansion and Q\\ appears as a subsidiary result which is inter- 
preted as the response of a plane detector parallel to a Lambertian uniformly 
distributed disc source. Circular apertures and sources are often considered 
in optics and radiation physics; and disc-shaped detectors are widely used in 
nuclear science (e.g. Neutron Activation Analysis). References to other works 
on this subject in the context of nuclear physics can be found in [6] and [7]. 
While the case of a point isotropic source has been treated to great extent 
[8,9,10,11,12,13], the case of a point cosine source has, to the best of our 
knowledge, attracted little attention, with the already mentioned exception of 
[3]. For these reasons, in the present work we extend the scope of the previous 
results by performing the calculation of the solid angle defined by a circular 
disc and a point cosine source pointing at an arbitrary direction, under the 
sole restriction that the disc lies in the half-space illuminated by the source. 

The solid angle subtended by a given surface at a point source located at the 
origin can be defined by 



n mrf = jj f(n)dn, (i) 

directions 
hitting surface 



where f(Q)dQ is the source angular distribution. In the case of a point cosine 
the distribution is defined in relation to some direction axis specified by the 
unit vector k and it is given by f(fl) = (f2 • k+ |f2 • k|)/(27r), the (27r) _1 
factor ensuring that < Q SU rf < 1- For A* — f2 • k it follows that f(fl) = 
{fi/n (fi > 0);0 (fx < 0)} so that the source only emits into the hemisphere 
around k. Because of this, the integration limits in the RHS of eq. 1 are to be 
determined from the conditions that /i > and that each included fl direction 
hits the surface. In the following we shall assume that the position of the disc 
is always chosen is such a manner that /i > for each point on the disc. This 
restriction greatly simplifies the calculation without, we believe, reducing the 
practical interest of the expressions. This is so because, in actual situations, 
the source is distributed over some planar surface; the source axis is coincident 
with the normal to the surface and then the restriction simply requires that 
the detector be held somewhere on the illuminated side of the radiating surface 
or eventually directly on the surface but does not intersect it. 



2 



2 Solid Angle Calculation 



To proceed it is advantageous to consider two coordinate systems (S' and 
S) with a common origin also coincident with the source position. In the S' 
system (fig. 1) the z' axis is aligned with the source direction k; the position 
of the disc center (C) is specified by h and R; and the symmetry axis of the 
disc is specified by the angles (5 and 7, (5 being the angle between the source 
and the disc axes. When working in the S system (fig. 2), the z axis is parallel 
to the disc axis, C is located by means of L and d; and k is given by angles (3 
and a. Due to the symmetry of the source it is possible to choose the x, y axes 
so that the y coordinate of C is zero in each coordinate system. Also from the 
symmetry of the source it is clear that the solid angle is an even function of 
a or 7 and in the following calculations we will thus assume that < a < it 
and < 7 < 7r. Finally, it is sufficient to consider < (3 < ir/2. 




Fig. 1. Notation used in the S' coordinate system. 




Fig. 2. Notation used in the S coordinate system. 
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The restriction that /x > over the whole disc is easily expressed in S' by 

h - r ■ sin /3 > , (2) 

r > being the disc radius. This is automatically satisfied, provided h > r, 
but otherwise, it further reduces the range of variation of [3 to 



< (3 < arcsin(/i/r) ,h<r 



(3) 



It is readily found that the coordinates of a point in S' and S (P' and P) are 
related by an orthogonal matrix (M): 



P' = M P , 



(4) 



where M is given by 



M( 7 ,M 



P = M T P' , 



(5) 



sin 7 sin a+ cos 7 cos (3 cos a sin 7 cos a— cos 7 cos (3 sin a —cos 7 sin (3 

cos 7 sin a— sin 7 cos /? cos a cos 7 cos a+ sin 7 cos /? sin a sin 7 sin /? 

sin /? cos a — sin /3 sin a cos (3 

(6) 



and its inverse is given by the transpose 



M T ( 7 ,/3,a) 



sin 7 sin a+ cos 7 cos /3 cos a cos 7 sin a— sin 7 cos /? cos a sin /? cos a 

sin 7 cos a— cos 7 cos /3 sin a cos 7 cos a + sin 7 cos /3 sin a — sin /3 sin a 

— cos 7 sin (3 sin 7 sin [3 cos /? 

(7) 



Let (x c , y c , z c ) and (x' c , y' c , z' c ) denote the coordinates of the disk center in S 
and S', respectively. As previously said, it is possible to choose the x, y axes 
in each coordinate system so that y c = y' c = 0. Starting in the S system, 
setting (x c , y c , z c ) = (d, 0, L) and fixing a and (3, one then uses eq. 4 to obtain 
(x' c , y' c , z' c ). The value of 7 is determined by imposing that y' c = 0, which gives 



tan 7 = {d sin a) / (d cos (3 cos a — L sin (3) , 



(8) 



or, for < 7 < 7r, 



4 



7 = i(/3,L,d,a) 



= arccos[(d cos (3 cos a — L sin (3)/\J (d sin a) 2 + (d cos (3 cos a — L sin f3) 2 ] . 

(9) 

The values of h and i? are obtained from 



h = z' c = h(/3, L, d,a) — d sin (3 cos a + L cos (3 , (10) 

and 



R = x 'c = y x 'c+y'c 
= R((3, L, d, a) 

= \J d 2 (sin 2 a + cos 2 a cos 2 (3) + L 2 sin 2 (3 — dL sin(2/3) cos a . (11) 
Eq. 9 can then be written as 

7 = arccos[(<icos / 9cosa — Lsm(3)/R] . (12) 

Conversely, starting in the S' system with (x' c ,y' c , z' c ) = (R,0,h), using eq. 5 
to obtain (x c ,y c ,z c ) and imposing that y c = 0, there results 



a = a(/3, h, R, 7) 

= arccos[(i? cos f3 cos 7 + h sin (3)/ ^J(Rsm'y) 2 + (i?cos/9cos7 + hsm[3) 2 ] . 

(13) 



L = z c = L(f3, h, R,j) — h cos (3 — Rsm (3 COS7 , (14) 

and 



d = x c = •y/x 2 + y 2 
= d ((3,h,R, 1 ) 

= ^(sin 2 7 + cos 2 7 cos 2 (3) + h 2 sin 2 (3 + Rh sin(2/3) cos 7 . (15) 
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Again, eq. 13 can be cast as 



a = arccos[(i?cos/9cos7 + hsmfi)/d\ . (16) 

The substitution of the expression for h (eq. 10) on the RHS of eq. 2 and a 
bit of algebra gives the restriction /i > expressed in terms of S variables: 

L > (r — dcosa) ■ tan (3 . (17) 
If (5 — 7r/2, the previous eq. can be rewritten as 

cosa >r/d. (18) 

The solid angle (eq. 1) is best calculated in S. With the notation described in 
fig. 2, it is seen that k = (sin (5 cos a, — sin (3 sin a, cos f3) and 



k • ft = cos (3 cos 9 + sin f3 sin 9 cos(a; + ip) , (19) 

where, for a given direction, 9 is the polar angle from the z axis and <p is the 
azimuthal angle in the xy plane measured from the x axis. The solid angle is 
then given by 

Because the position of the disc is such that // > (eq. 17 or, equivalently, 
eq. 2), the integrations limits are determined only by the condition that each 
included (9, </?) direction hits the detector. 



d = OC 
P_= OA 




Fig. 3. Integrations limits for <p and definitions of p±. 
Referring to figs. 2 and 3, it follows that y? max = — Win = v^o = arcsin(r/c?) , 

#max = arccos(L/y / L 2 + p^_) , 
6» min = arccos(L/ v /L 2 + , 
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where 

p± = d cos ip ± \Jr 2 — (dsinip) 2 . 
Then, 



r + lpO rfmax 

f2 = 7r~ (cos /3 / / cos 6 sin 6> d6dip+ 



mm 

"Vax 



where 



r~r" i 

sin (3 / cos(a + <p) sin 2 6» dOdip) , (21) 

•> —<fO Jdmin 

= cos /3 + sin f3 S~2j_ , (22) 



(- + 1,00 fffmax 

Q^ti- 1 / cos9 sm9 dOdip (23) 

J-(po ./0 min 



and 



fi ± = vr" 1 / cos(a + (p) sin 2 dedip. (24) 

J-(pO "'Smin 

The values of Qy and f2j_were previously obtained. From [3, eq. 29], or [2, eq. 
33], 



1 1 (r 2 — J 2 — T 2 \ 

n, (M , rH - [1 + -_ L__J 1 (25) 

and from [1, eqs. 20 and 55], 



n ± (L, d,r, a) = cos a-^[ 1 - 1] , (26) 

2a VI - 



V . 1 1 '~ 

where 



m = 2rd/(L 2 + d 2 + r 2 ) . (27) 

In eq. 26 the absolute value of L is used since L is restricted only by eq. 17 
and can thus be negative. 

2.1 Special cases 



We emphasize two cases: (i) the axis of symmetry of the disc is parallel to that 
of the source ((3 = 0); and (ii) the center of the disc is located on the source 
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axis (R — 0). 



(i) When (5 — 0, from eqs. 9, 10 and 11 there results that 7 = a, h = L and 
R = d. The same results could of course be obtained from 13, 14 and 15. 
The restriction imposed by either of the equivalent eqs. 17 or 2 gives L > 0. 
Using eqs. 22, 25 and 27 it is seen that the solid angle is independent from 
a as expected and 

n(P = 0)=n\\(L,d,r) (28) 
It was shown in [2] that Q\\ is continuous except for L — > since 

fi„(L -> 0) = {0 (d > r), 1/2 (d = r), 1 (d < r)} . (29) 

(ii) When R = 0, eqs. 13, 14 and 15 give a = , V7, L = h cos (3 and d = h sin /3. 

Using eqs. 22, 25, 26 and 27, a little algebra yields 



= 0) 



cos 3 r 2 



yV + h 2 + 2rh sin J) (r 2 + /i 2 - 2r/i sin J) ' 

Q(R = 0) is then independent from 7, as expected. 
Eq. 30 can be approximated by 



(30) 



n(R = 0) = (£) 2 cos/5 



-J cos 2/5 + 



,r«/i. (31) 



In the case where h — r, it is straightforward to simplify eq. 30 to: 

n(R = 0,h = r) = 1/2 ,V/3 . (32) 



3 Results 



Here we present sample plots of the solid angle. In the following we choose to 
work with S' parameters (3, h, R, 7) and consider throughout a disc of radius 
r — 1. We first address the case /? = 0, so that Q = Q\\ , L = h and d = R. 
As said before (section 2.1) f2 is not continuous as L — > 0. This is illustrated 
in fig. 4 where f2 is plotted as a function of h(L) for different i?(d) values. 

In fig. 5 the situation where the disc center is on the source axis (R = 0) is 
represented. Here, Q does not depend on 7 (section 2.1) and plots of the solid 
angle as a function of (5 are shown for different h values , including the special 
case h — r, where Q = 1/2, regardless of (5. To avoid restricting the range 
of variation of 3 (see eq. 3), the values of h were chosen such that h > r. In 
figs. 6, 7 and 8 the effect of offsetting the position of the disc center from the 
source axis is shown for three offset values (R — h/A, R — h/2 and R = 3 /Ah) 
and for two values of 3 (vr/12, n/Q). As argued before, the solid angle is an 
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Fig. 4. Solid angle subtended by a disc of radius 1 with symmetry axis parallel to 
that of the source. 
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Fig. 5. Solid angle subtended by a disc of radius 1 and center located on the source 
axis. 

even function of 7 and it is seen that for the larger tilting angle (f3 = 7r/6) the 
dependence on 7 is enhanced. 
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Fig. 6. Solid angle subtended by a disc of radius 1 for fixed h = 2. 




Fig. 7. Solid angle subtended by a disc of radius 1 for fixed h = 3. 

When plotting Q as a function of 7 while holding all other parameters constant, 
there can be a zero, i.e. a value of 7 such that fi(7o) — 0- Looking at fig. 2 it 
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Fig. 8. Solid angle subtended by a disc of radius 1 for fixed h = 5. 

is clear that f2 = if L = and d > r. Setting L = in eq. 14 and solving 
for 7 gives an equation for the zero, 



70 = arccos[/i/(i?tan/3)] , (33) 

which has no solution for h/(R tan (3) > 1. Thus, this qualitative feature is de- 
pendent on the sign of h — R-tanfl, which actually distinguishes the situations 
where the disc always presents the same face to the source from those where 
the face presented depends on the value of 7. This is schematically explained 
in fig. 9 where the effect of changing 7 from to tt is shown as seen in S' (see 
also fig.l), when looking along the yf axis. If h — R ■ tan/3 = (fig. 9a), the 
source always 'sees' the lower face of the disc and for 7 = one has L = and, 
consequently, Q — 0. If h — R • tan (3 > (fig. 9b) the source always looks at 
the lower face of the disc but f2 is never zero. Finally, for h — R • tan (3 < (fig. 
9c), it is seen that, as the disc swirls with increasing 7, the upper face is first 
shown (e.g. 7 = 0) and then the lower face is presented (e.g. 7 = tt), which 
means that Q = at some point in between. This behaviour is illustrated in 
fig. 10, for h — 2 and (3 = tt/3. The zero shows up for R> h/ tan/3 ~ 1.16. 

One should notice that in the preceding discussion it was implicitly assumed 
that d > r when 7 = 7o- We now proceed to show that eq. 2, in the strict form 

h - r ■ sin f3 > , (34) 
guarantees that d > r. Using eq. 33 to eliminate h in eqs. 34 and 15 gives 

cos 7 • Rj cos (3 > r (35) 

and 



d = Rj cos (3 ■ y cos 2 7 + sin 2 7 cos 2 (3 . (36) 
respectively. Since J cos 2 7 + sin 2 7 cos 2 (3 > cos 7 therefore d > r. 
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(a) (b) (c) 

Fig. 9. The effect of changing 7 as seen from a side view in the S' system for: (a) 
h = R- tan((3) , (b) h > R ■ tan(0) and (c) h< R ■ tan{/3). 
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Fig. 10. Plots of Q as a function of 7 showing a zero for h < R ■ tan((5). 



4 Summary and Outlook 



Analytical expressions for the solid angle subtended by a circular disc at a 
point cosine source were obtained , under the single restriction that the disc 
is located in the half-space illuminated by the source (eq. 2 or eq. 17). It was 
shown (eq. 22) that the solid angle can be decomposed into the combination 
of two components and Qj_) corresponding to the situations where the 
symmetry axis of the disc is parallel (eq. 25) or orthogonal (eq. 26) to the 
source direction. Q can be calculated relatively to two alternative coordinate 
systems (S' and S) shown in figs. 1 and 2. The parameters pertaining to each 
system ({7, h, R, r)}, {a, L, d, ((3,r)}) are related through eqs. 9, 10 , 11, 
13, 14 and 15. 

A similar calculation to that presented here can be performed for the solid 
angle defined by a cylindrical detector. A work where we report this latter 
result is in preparation. 
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